Growth of graph states in quantum networks 
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We propose a scheme to distribute graph states over quantum networks in the presence of noise 
in the channels and in the operations. The protocol can be implemented efficiently for large graph 
sates of arbitrary (complex) topology. We benchmark our scheme with two protocols where each 
connected component is prepared in a node belonging to the component and subsequently distributed 
via quantum repeaters to the remaining connected nodes. We show that the fidelity of the generated 
graphs can be written as the partition function of a classical Ising-type Hamiltonian. We give exact 
expressions of the fidelity of the linear cluster and results for its decay rate in random graphs with 
arbitrary (uncorrelated) degree distributions. 
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I. INTRODUCTION 

Quantum networks [3, where nodes with (lim- 
ited) quantum storage and processing power are coupled 
through quantum channels, are becoming a focus of in- 
terest in quantum information. The first motivation is to 
extend the paradigmatic bipartite quantum communica- 
tion applications, e.g., quantum teleportation [3| or more 
prominently quantum key distribution [1] , to a multipar- 
tite setting, where such bipartite protocols can be accom- 
plished between arbitrary nodes of the network. The pos- 
sibility of networks with non-trivial topologies can give 
rise to new phenomena and to applications that exploit 
the multipartite correlations. Of course, this puts for- 
ward a big variety of theoretical and technological chal- 
lenges which can be addressed in short term. Indeed 
the first steps have already been taken towards a quan- 
tum network in a first realization with two distant nodes 
that can store and interchange quantum information in 
an efficient and reversible way Q. Lastly, such physi- 
cal realizations are in principle scalable and hence open 
the door to perform highly controllable experiments on 
many-body phenomena, study multipartite entanglement 
Q , and could eventually perform more complex tasks like 
some implementation of the quantum Google page rank 
or general distributed quantum computations. 

A central task in quantum networks is to devise strate- 
gies to distribute entanglement among its nodes. Linear 
networks have shown to be useful for long-distance bi- 
partite entanglement distribution by means of quantum 
repeaters [|[ (see also its measurement-based implemen- 
tation [9j]). However, the study of entanglement distri- 
bution over higher dimensional networks is in its infancy. 
Current results show that in some scenarios the network 
topology can b ring interesting effects like entanglement 
percolation [1CH16I ] that lead to new approaches to the 
problem. This effect shows up in non-ideal (but still not 
fully realistic) scenarios or for some particular network 
geometries. 

In this paper we address the problem of distributing 
graph states in the realistic scenario of noisy network 



channels and a small, but non- negligible, amount of noise 
in the local operations. Graph states [ij} are a large 
family of multipartite entangled states that, although 
they can be efficiently described with relatively few pa- 
rameters, have a rich variety of features. For instance, 
they include paradigmatic states like GHZ, cluster states, 
codewords of error-correcting codes (Lsj . they provide 
some novel quantum communication applications, like 
secret entangled-state distribution [l^, they are use- 
ful in the study of non- locality [2fJ |. and most impor- 
tantly they include states which are universal resources 
for measurement-based quantum computing [2l|, [H| . 

Our goal here is to propose a protocol to create a 
large graph using an underlying network of noisy chan- 
nels. The protocol should tolerate channels and oper- 
ations with errors, scale efficiently with the size of the 
network, and work for any network topology, and inpar- 
ticular for complex networks. Complex networks [23l425| 
underlie many natural [26|, [27j], social [28l . [29j and ar- 
tificial systems [13, HH where different parties interact. 
Their non-trivial structure is the source of features, like 
the existence and shape of an important fraction of highly 
connected nodes [HJ and the tendency of nodes to clus- 
ter together (28|, that are often observed in real graphs 
and have a deep impact in their performance. Complex 
networks are of particular importance in communication 
infrastructures, as most present telecommunication net- 
works like Internet have a complex structure [13, [HI- 
However, the rich and intriguing properties of these type 
of networks are still quite unexplored in the quantum 
setting (see nonetheless [H El HI, H ) • 

Channels linking separate nodes in a network are typ- 
ically noisy, and pose the main caveat to the creation 
of distributed multipartite entangled states with hig h fi- 
delity. To overcome this, there exist bipartite [35ll36l| and 
multipartite |37H33 | entanglement purification protocols 
that allow either to generate highly purified Bell pairs, 
which can latter be used to telcport an arbitrary graph 
state, or to directly purify the desired graph state. These 
recursive protocols tolerate a reasonable amount of noise 
in local operations, but require a number of initial copies 
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that grows exponentially with the size of the state. Other 
proposals do not use postselection, making the purifica- 
tion efficient in terms of the size of the graph state, but 
come at the expense of a stricter noise threshold [40| . 
In Ref. [4lj |. an entanglement pumping without postse- 
lection is also used to obtain efficient purification when 
constructing the graph state edge by edge. In different 
approaches, the graph state is created by a probabilis- 
tic^grow th using non-deterministic entangling operations 

Here we investigate the advantages of generating and 
purifying small GHZs that reproduce the local structure 
of the network, and merge them in order to distribute a 
network-wide graph state. Since GHZ have a fixed size 
which depends on the degree of each node and is thus 
independent of the size of the network, this protocol is 
efficient in the size while it still maintains the high thresh- 
olds of the recurrence schemes. We benchmark this pro- 
tocol with two other protocols that generate high-fidelity 
bipartite states between a node and the rest of the net- 
work, which are then used to distribute a locally gener- 
ated graph state. We use the fidelity of the graph state 
as a figure of merit to compare the three protocols. The 
fidelity decays exponentially with the size of the network 
for a constant level of noise, and hence we also use its 
decay rate. It turns out that both quantities can be un- 
derstood as the partition function and free energy of a 
thermodynamic system, respectively, and thus standard 
methods of statistical mechanics are readily used. 

The paper is structured as follows. First, in Section HT1 
we introduce graph states and the noise model we con- 
sider. In Section HTT1 we present the three protocols and 
the two figures of merit — the fidelity of the graph state 
and its decay rate — , and relate them with an analog 
partition function and free energy of a thermodynamical 
system. Then, in Section IIVI we apply the protocols to 
the creation of a linear cluster state, for which we ob- 
tain exact results, and of a graph state associated to a 
complex network. We conclude in Section IV1 



II. DEFINITIONS 



A graph state [l7( is a quantum state associated to a 
graph G, where vertices correspond to qubits and edges 
to interactions. A common description of graph states 
is the stabilizer formalism. For a graph G, there are 
N = \V(G)\ stabilizer operators (one for each vertex it) 
defined as 



K u = X u 



n z v . 



(1) 



Here X, Y and Z are the Pauli operators and the 
subindcx denotes on which qubit they operate. A pure 
graph state \fi) G , with fi — {0, 1}^, is a common eigen- 
state of all stabilizer operators with 

k% |/x) G = (-ly- \ti) G v^y, 

where ji u is the u-th component of fi. The set {|m) g } 
form the graph state basis. We will omit G when the 
graph is clear by context. 

An alternative description is by means of the interac- 
tion picture, where |0) G is created by preparing all qubits 
in the state |+) and then applying a CPHASE Mj „ for every 
edge (u, v) € E, 



|0> G = 



Yl CPHASE^ |+) 

(u,v)£E 



Since Z u anticommutes with K u and commutes with the 
rest of K v , v =^ it, any graph-basis element can be ex- 
pressed as 



G ■ 



A. Graph states 



Throughout this paper we consider mixed graph states 
diagonal in the graph state basis, pa — |m) {/-*!• 

Any mixed graph state can be brought to this form by 
local depolarization (38[. We will work however in the 
stabilizer basis. Let us define K x = Y\ u£V (K v ) v , where 
x = (xiX2 ■ ■ -Xn) € {0, 1}^. Since {K x } form a com- 
plete set of commuting observables, a diagonal graph 
state can be expressed as pa — J2 X where 
(K x ) = £ p A M (-l) M a5 . The effect of the Pauli operators 
on K, r is 



A graph G — {V(G),E(G)} is an ordered pair of sets: 
the set V(G) of vertices, or nodes, and the set E(G) 
of edges, or links, whose elements are (unordered) pairs 
of vertices and represent the connections between them. 
The neighborhood of a vertex it e V(G) is the set of ver- 
tices connected to it, M U {G) — {v : (u, v) £ E{G)}, and 
the degree of that vertex is the number of its neighbors, 
k = \Af u (G)\. If the graph is directed, then the elements of 
E(G) are ordered, and (u, v) 6 E(G) is a directed edge 
from it to ti. In this case, the incoming and outgoing 

neighborhoods of it are Nu n \G) = {v : (v,u) 6 E(G)} 
and A^° ut) (G) = {v : (u,v) £ E(G)}, respectively. 



Z U K X Z U = (-1) X *K X , (2) 
X U K X X U = \{ (-l)**K a , (3) 

Y U K X Y U = {-l) x " 11 (-1) X *K X , (4) 

so this unitaries map diagonal graph states into diagonal 
graph states. A CPHASE^ adds an edge between it and 
v, if there where not connected, or removes it, if the edge 
already existed. The effect on K x is 

cphase u ,„.R: x cphase u ,„ = K, (Z v ) Xu (Z u ) x " . (5) 
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The action of Pauli measurements can also be easily 
described in this formalism as a transformation of the 
graph (up to some local unitaries). Measurement of Z 
simply disconnects the measured qubit from the rest of 
the graph, while X and Y transform the neighborhood 
of the measured qubit and then disconnect it. In terms 
of the stabilizer operators, the measurement of Z u com- 
mutes with all K v , for v u, and anticommutes with 
K u . Thus, 

[I+(-l) m Z u ]K x [I+(-l) m Z u ] 

= [I+(-l) m Z u ]K x S ^, (6) 

where m = {0, 1} labels the measurement outcome 
{+1,-1} respectively. After tracing out qubit u, the 
new stabilizer is 

(-l)"^e^^K x So, x ^ (7) 

where the new K' x corresponds to a new graph G' ob- 
tained form G by removing vertex u and its attached 
edges. The spurious phase factor in (0 can be cancelled 
by applying a unitary (IlueJV %v) ■ Similarly, measure- 
ments of X u or Y u also result in the disconnection of the 
measured qubit, but in these cases the remaining graph is 
transformed by local complementations of the neighbor- 
hood of u, as described in [13]. We will be more explicit in 
the concrete cases where we apply these measurements. 

Our figure of merit will be given in terms of the fidelity, 
which for a graph-diagonal state can be written as, 

F% = G {0\ Pg \0) g = A = ^ ]T ( K S) • (8) 



B. Network and noise model 

We consider a network where nodes are spatially sepa- 
rated. Some of the nodes are connected by a link, which 
we model as a noisy depolarizing channel on one qubit 
with error parameter p c , 

TM=(l-p c )[/] + |^[^]. (9) 

i=Q 

Here ai are /, X, Y, Z for i = 0, 1,2,3, respectively. A 
Pauli measurement on u is modeled as a perfect measure- 
ment preceded by a depolarizing channel Tj with error 
probability p\ on that qubit, 

3f tt, = (l-Pi)[i] + | X>i U) ]. (10) 

t=0 

A noisy two-qubit gate (e.g., a CPHASE or a CNOt) on 
qubits u and v is modeled as an ideal gate followed by 
the two-qubit depolarizing channel on u, v with error pa- 
rameter p 2 , 

i,j=0 



Gates can only be applied locally, i.e., on qubits within 
the same node. 

In the stabilizer basis, the effect of each noise source 
is easily tracked: it multiplies each stabilizer element 
by a coefficient, and keeps the graph state in diagonal 
form. A noise T 2 changes the sign of the coefficient 
accompanying K x with probability P2/2 if it acts non- 
trivially on u or v, i.e., if exists at least one x a — 1 for 
sGiiUd L)Af u UAf v . Thus, each T^"' 1 "' adds a coefficient 

(1 — p 2 ) e( - x "' x ^' x ^^ x ^\ (12) 

where 6(x) = 1 — 5(0, x) and xj^ u = (x Vl . . . x Vk ) for all 

v a G M u . Noise t[ u ^ behaves similarly. In this case, the 
sign is changed with probability pi/2 unless x u = and 
@ a x a = for a G J\f u , so the multiplying coefficient is 

(l_p 1 )*(*«.©. e Af««-). (13) 

Note thus that these noise sources affect the qubits on 
which the gates act, plus their neighbors. 

Finally let us point that we do not associate any noise 
to the local "correcting" unitaries performed in order to 
bring the post-measurement states to a standard graph 
form. We refer to this unitaries in our protocols to sim- 
plify bookkeeping, but their action can be pushed forward 
(or commuted) till the end of the protocol, and hence the 
resulting state is exactly equivalent as a resource of en- 
tanglement. 

III. PROTOCOLS 

We propose a protocol to distribute a graph state with 
the structure of a general quantum communication net- 
work, of arbitrary topology, associated to graph G. Noise 
p c in the communication channels is considered to be rel- 
atively high, so some sort of purification or error correc- 
tion is in order. The protocol uses the structure of the 
network to distribute several copies of small subgraphs 
between neighbors, and purifies them by means of mul- 
tipartite purification. These subgraphs are GHZ states, 
which are associated to a star graph with a central node 
of degree j, connected to j leaves of degree 1. In order 
to benchmark our protocol, we also consider two refer- 
ence protocols that distribute bipartite states between 
a central node, which locally creates the desired graph 
state, and the rest of the network. The bipartite states 
are then used to teleport the locally created graph state. 
Since this central node may not be directly connected to 
the rest of the network, quantum repeaters Q are used 
to establish purified bipartite states between this node 
and all the network's nodes. 

In all cases, we consider the same multipartite purifi- 
cation protocol for bicolorable graph states described in 
Ref. [38( • The noise threshold of this purification scheme 
depends on the degree of the central qubit of the state 
that is purified, and goes from w 0.06 for a bipartite state 
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Figure 1: Bipartite A protocol. The upper green area rep- 
resents the "central" node, where the local graph state with 
the structure of the network is generated. The lower green 
circles are nodes in the network. Small, gray dots correspond 
to qubits, and lines connect neighbors. Rectangles in red in- 
dicate Bell measurements involved in the teleportation of the 
local graph states, which are implemented as CPHASEs be- 
tween u and u a and X measurements at u, u a . 



(i.e., a GHZ with central degree 1) to « 0.02 for a GHZ 
of central degree 9. This poses a limit in the maximum 
degree of a network for which the protocols can be used. 
Hence, we work under the assumption that the errors in 
the local operations are low — as compared to the (finite) 
degree of the subgraphs, jp <C 1 — and that the channel 
noise can be significant, but also low enough so that the 
minimum threshold fidelity required for the purification 
protocol to succeed can be attained. At first order in 
Pi = p2 = p, the output state fixed point of this purifi- 
cation scheme is given by (see Appendix [XJ 




1-p 



+ X a (j + 1) 



(14) 



where a is the central node and b the leaves of the GHZ, 
j is the number of leaves and the overline in x a represents 
the bit-complement, = x a © 1. 

Note that the extensive use of quantum repeaters in 
the first two protocols renders them extremely inefficient. 
Nevertheless, we will find that the performance in terms 
of attainable fidelity is still comparable to that of our 
more efficient subgraph protocol. 



A. Bipartite A protocol 

In the first protocol, purified, but still noisy, entan- 
gled states are created between a central node and all 
the others by means of quantum repeaters. This central 
node then teleports a locally generated state to all the 
other nodes. For simplicity, in protocols A and B we as- 
sume that the central node is sending all qubits through 



purified channels, including his own. This adds a source 
of noise that would not be strictly necessary, but its effect 
is small for large enough networks. In cases, like complex 
networks, where the network is not necessarily connected, 
it is understood that there is a central node for each con- 
nected component distributing the corresponding graph 
state. 

The local state pa, which mimics the structure of the 
network, is created by initializing N qubits in the state 
stabilized by X and then applying CPHASEs between 
neighbors. Then, the noisy Bell states between the cen- 
tral node and all the others are used to teleport the cor- 
responding qubits (see Fig. [1]). We label qubits in the 
local graph as u £ V(G). For each node in the network, 
there is a bipartite state p 9u of two qubits, u a , Ub 6 V(g u ) 
and one edge (u a ,Ub) € E(g u ). Qubit u a belongs to the 
central node, while Ub is in the corresponding node in the 
network. In order to account for the errors in the telepor- 
tation Bell-measurement, we implement it by a CPHASE 
on (u,u a ) followed by X measurements on u and u a . 

Without taking noise into account, the state before 
teleportation is 



with 



PG 



PG® p 3n 

ueV(G) 



and pg u 



(15) 



The graphs g u correspond to the bipartite states with 
qubits u a and Ub used for teleportation, i.e., — 
X Ua Z Ub . The action of the CPHASE„ jUo affects only the 
stabilizers K® -)• K%Z Ua and Kjfc K^Z U . The mea- 
surement of X u anticommutes with all K^f , v € Af u (G), 
and with K^ U Z U , while that of X Ua anticommutes only 
with K^Z Ua and K%?. Thus, each term changes to 



[1 + (-l) m ^X u J [1 + (-ir«Aj K a 



(16) 



where m u and m Ua are the 
comes. Tracing out qubits u 



measurement out- 
and u a and cor- 
recting the state with (Z Ub )" Lu (X Ub ) " a we obtain 
K ^o,x Ua <s® ve ^ u xJa,x u ®x Ui , with K' x associated to a 
new graph G' where vertex u has been substituted by 
vertex Ub- Teleportation of all local qubits results in the 
desired distributed state. 

Noise can now be introduced as the multiplicative fac- 
tors in front of the stabilizer elements of (fT"5|) . Here, 
the order of the CPHASE gates used to generate the lo- 
cal graph matters, as qubits receive noise from gates 
performed at neighbors which are already connected to 
them. There is thus a noise (1 — p\) Xu correspond- 
ing to the preparation of each node in X Ul and a 
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Figure 2: Bipartite B protocol. The upper green area rep- 
resents the "central" node, where the local graph state with 
the structure of the network is generated. The lower green 
circles are nodes in the network. Small, gray dots correspond 
to qubits, and lines connect neighbors. Dashed lines, corre- 
sponding to the edges of the local graph, are created after 
solid lines, corresponding to the distributed bipartite states. 
Squares in red indicate X measurements used to teleport the 
local graph states. 




Figure 3: Subgraphs purification and connection protocol. 
Green circles are nodes in the network. Nodes on the right 
are the outgoing neighborhood of it, while those on the left 
are the incoming neighborhood. Small, gray dots correspond 
to qubits, and lines connect neighbors. Lines and squares in 
red indicate CPHASEs and measurements involved in the con- 
nection of subgraphs. 



(1 - p 2 ) e ( x «' x '» x J ! ?«' x J ! ?J for each edge in the local state, 
where the tilde in Jv u labels that we consider the neigh- 
borhood of u at the moment the CPHASE Uj „ is per- 
formed. Additionally, and using that x Ua = @ ve j^ u x v 
and x Ub — x u , in each teleportation the CPHASE in- 
troduces noise as (1 — p 2 ) et ^ Xu ' XN ^\ and the measure- 
ments as (1 — Pl y x ^ + ®v€Ar u x " _ Finally, there is a 

(K^)®«^ Xv (K^) x A term from the purified Bell 
states. This results in a final distributed state p — 

(k*)= n [(i-pi) ai "((^)®- 6v - B "W) B ") 

uev 

(l_p 2 )«(^.^«)(l_ pi )^+©»e^^ 

n (! -p 2 ) Hxu ' xv,a> *"' iB *» ) , (17) 

where the tilde over n states that the edges are intro- 
duced in a certain order. An explicit expression for the 
Bell state correlators for the postselection purification 
protocol at first order in p can be obtained from f| 14[) : 

{{K u J^{K Ub ) x ->) ~ 1 - (x Ua x Ub +2x Ua ) P . (18) 

B. Bipartite B protocol 

The previous strategy can be improved by directly con- 
necting the local ends of Bell pairs by means of CPHASEs, 
and then teleporting the local graph state performing 
only one X measurement per node — instead of the var- 
ious sources of noise induced by the Bell-measurement 



(cphase and two X measurements). The initial state is 
now ®«ev(G) P9v.i wner e qubits are labeled as u if they 
belong to the central node and Ub if they are in the dis- 
tributed nodes (see Fig. [2]). The sources of noise are: 
the purity of Bell pairs, CPHASEs used in the prepara- 
tion of the local state and the measurement involved in 
each teleportation. As in the previous case, the order of 
CPHASEs is important: each contributes to a noise with 
(1 - p 2 ) 6 ( x «' x '" x »b> x '>b' ,B tf*> a! *rJ, Note that the CPHASEs 
are performed on qubits that are already connected to 
nodes in the network and thus affect the distributed 
qubits Ub- The calculation is similar as before and re- 
sults in a final distributed state with correlators 

J^J (T — p2^®( Xl " Xv '®™eM u x ™-®weMv x w x Ar u > x jVv) _ 
(u,v)eE 

(19) 

C. Purify subgraph and merge 

In this strategy, subgraph states of small size (N inde- 
pendent) are distributed and purified, and then intercon- 
nected locally at each node to form the desired structure. 

The protocol follows two steps (see Fig. [3]). To each 
node of degree k we assign an outgoing neighborhood, 
with j nodes, and an incoming neighborhood, with i 
nodes (k = j + i). First, each node prepares a GHZ 
of size j + 1. A GHZ graph state with j + 1 qubits is 
associated to a star graph g u , i.e., a graph with a cen- 
tral node and j external nodes, called leaves. Then, each 
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qubit corresponding to a leaf of the GHZ is sent to one of 
the outgoing neighboring nodes through the depolarizing 
channels. Several copies of this distributed state are cre- 
ated and th en p urified using the bicolorable graph state 
protocol of [37|. The final purified subgraph state with 
central qubit u is 



pu=± e ((w n (*£)"-) 

(xtr n ( KS ui) Xuw > (2°) 

where the correlators of the fixed point of the purification 
scheme are given in Eq. (fl4|) . Here u w denotes the leaf 
qubit that has been sent to w £ 7V"i° ut ^(G). At the same 
time, each node receives i — k — j leaves corresponding 
to the GHZ states created at its incoming neighborhood. 

In the second step, each node connects the central 
qubits of their GHZ state with the leaves they have re- 
ceived. The connection is made by performing a CPHASE 
between the two qubits u and v u and a Y measurement of 
the received qubit v u . The effect of the Y measurement 
is to add an edge between the two central nodes u and v 
and destroy de measured qubit v u . 

The action on the stabilizers can be seen in a sin- 
gle connection example between a central qubit u of 
a subgraph and a leaf qubit v u of an incoming sub- 
graph (see Fig. [3]). We focus on the stabilizer elements 
(K§ v ) Xv (K^) Xv ^ (K^ u ) Xu , as all the others remain un- 
changed. The action of a CPHASE between u and v u is 
K%l -> Kg»Z u and ->• K$»Z Vn . The measurement 
of Y Vu anticommutes with the new stabilizer operators at 
v, v u , and u, so each term changes to 

[1 + (-ir%j (K^(K^z u r^(K^z Vu r- 

3o,x v ($Xv u ®x u - (21) 

Tracing out v u and correcting the state de- 
pending on the measurement outcome m with 
exp[(-l) m if Z v ] exp[(-l) m if Z u ], we obtain 



(K' v r<>(KZ«Z v r«6 



where K' v is the stabilizer of v with qubit v u changed to 
u. 

Noise added by the measurement enters as 
(1 — p{Y u ® Xv , where we used Eq. (fl~3| together with 
Noise by a CPHASE between one of the in- 

, 9(X U1 X V ,JE ( ou t) ~ (in) ) 

coming leaves v u and it adds (1— P2) Mu Uu , 

where Nu^ is the incoming neighborhood that has al- 
ready been connected to u. 

Finally, note that if a node has j = it does not need to 
prepare any GHZ state. In this case, one of the incoming 
leaves is used as the qubit to which all the other leaves are 
connected. If there are no incoming leaves, then the node 
is isolated (and is thus prepared in the state stabilized 



by X). This means that the contribution of the GHZ 
distribution and purification at a node u is 

wr- n (kziy^A if j > 0, 

1 if j = 0, i > 0, 

(1- Pl ) x " ]£j = i = o. 

All in all, the final correlators for the subgraph protocol 
can be written as, 



(k x ) = n \ (Kg-) n [a-pi) 3 



(22) 



Here, denotes the incoming neighborhood M™ if 
j > 0, or the incoming neighborhood minus the first in- 
coming neighbor, Af™\vi, in the case j = where the 
first neighbor v\ is used as the qubit to which all remain- 
ing edges are connected. 



D. Fidelity 

The fidelity for a distributed graph G is F G — 
W Sec (^)j and is in general hard to calculate. Setting 
Pi = P2 = P, all errors from the CPHASEs and measure- 
ments enter as factors 

(l_p)h(«0 =e -W«0 (23) 

where we have defined ft = — log(l — p). To first order 
in p, the correlators of the purified graph states ([M)) can 
also be written in the same form, 

1 - h(x)p ~ e~ ph ^ ~ e - Wx) . (24) 

Hence, the correlators of the generated graph state can 
be written as (K x ) = exp[— (3H G (x)}, where H G (x) is 
the sum of the different noise terms h(x). The fidelity 
thus resembles the partition function of a system with 
Hamiltonian H G (x) and inverse temperature /?: 



-,7V X/ 



(25) 



The "Hamiltonian" H G (x) can be expressed as the 
sum of many-body n-local Hamiltonians of the form 



%2 j • * • j 


) = 1- 




£*=1 X b 




2 





Xi X2 ■ ■ ■ X n 

j j N 



6=1 



(26) 
(27) 

(28) 



6=1 



The last term ©£=1 Xb m ^1- is a j-body interac- 
tion term. By recursively using Eq. (|27p . this term can be 
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seen to be equal to the sum over all index permutations 
of J2 a odd x i ■ ■ ■ x a x a +i ■ ■ ■ xj. That is, we have rephrased 
our problem of computing the fidelity of a distributed 
large graph state as that of computing the thermal prop- 
erties of a classical many-body Ising-type system, where 
the indices x take the role of classic spins. The corre- 
sponding Hamiltonian will inherit the topology of the 
underlying graph and its precise expression will depend 
on the graph-growth protocol used. 

We are interested in the rate at which F$ decays, 



(29) 



where Nf§ is the analog of the free energy of the sys- 
tem. A good reason to study this quantity is that in 
statistical systems such as a complex network (which is 
modeled as an ensemble of graphs, each associated to a 
probability — see Section HV Bp . the partition function (fi- 
delity) in itself is not an extensive quantity, while the free 
energy is typically extensive and self-averaging (see, for 
example, Ref. P- 188). 

We can further exploit the statistical physics analogy 
and apply the known methods and understanding to com- 
pute the rate at which the fidelity decays for the differ- 
ent proposed protocols. We are interested in a regime 
where the noise in the operations is low, which corre- 
sponds to the high temperature limit. In addition, in the 
cases under study, each spin (or node) is effectively cou- 
pled to several spins, either as nearest- or second-nearest 
neighbors. These are conditions for which mean-field ap- 
proximation is very well suited: Eqs. ([26)) to (l28|) can 
be linearized using the standard mean-field approxima- 
tion to express the Hamiltonian as the sum of a constant 
term plus linear terms in x u . For this purpose we take 
x u — > s+ 5 U (and x~^ — > 1 — s — 5 U ) where s is the value of 
the mean-field and 6 U are the arguably small fluctuations 
of x u around its mean value. Keeping only the linear 
terms in the fluctuations, the different terms present in 
the Hamiltonian become, 



\X± , X2 i • • • , X 



„) = i-a- sr + a-s)"- 1 ^ 



(30) 

x u © x v = 2s(l - s) + (1 - 2s)(5 u + S v ) (31) 



1 - s 



6=1 



1 -S ( 1 



{^[1- (1-2^] 



+(i-2.sr i J> 



6=1 



^[l-(l-2^] 



(32) 



With this linearization, the new mean-field Hamilto- 
nian takes the form H MF = J2 U A U + Y^ u B u x u , where 



A u and B u are functions of s. Hence, the sum over x, 
i.e., "configurations", in the fidelity can now be carried 
out with ease, 



-PA Ue -PB/2 CQsh 



II' 



and the decay rate becomes (3f Mb , with 



(33) 



(34) 

The value of s is found by adding an artificial linear term 
to the Hamiltonian (playing the role of an external mag- 
netic field), changing B u — > B u + £, and requiring con- 
sistency in the definition 



Of 



MF 



One is hence left with the trascendental equation for s, 

1 1 ^ BB U 
— — > tanh — — . 
2N ^ 2 



(35) 



Its solution s* can be substituted back in the expression 
for / MF ([3^| to obtain the desired result. In the cases 
of interest here, B u is some polynomial which remains 
bounded for all values of s € [0,1]. Hence, to leading 
order in p we can approximate tanh ^rp- , arriving 

to 



1 1 

~ 2 -- 2 vB u \ 



(36) 



for small enough p. 



IV. NETWORK EXAMPLES 
A. Closed linear cluster 

We first study the case of a one-dimensional network, 
in which a linear cluster state is created. This case is 
remarkable because we can compute the exact fidelity 
for any cluster size. For symmetry, we consider a closed 
linear cluster state, where all nodes have degree 2. The 
order of the CPHASEs in the creation of the local graph in 
protocols Bipartite A and B, and the size of subgraphs in 
protocol Subgraphs, gives different results for the fidelity. 
For simplicity, in Bipartite A and B, we consider CPHASEs 
applied to successive nodes. Except for the noise of the 
first gate (which affects only two nodes) and of the last 
one (which affects four nodes), all the gates contribute to 



8 



Figure 4: Subgraph protocol in a linear cluster. Top: SI, 
using subgraphs of degree 1. Bottom: S2, using subgraphs of 
degree 2. 



the noise of three nodes. Each of the correlators is thus 

(k x ) ~ n (i 

(1 — p 2 ) 29 ^ x "~ 1,:c '" Xu+1 ^ 

(TO*- 1 ^ (37) 



for the Bipartite A protocol and 



{K x )c± JJ(l- Pl )^(l- F2 



(x u _i,a;„,x„+i,x„+2) 



(38) 



for the Bipartite B protocol. Recall that the correlator 
of the purified subgraphs is given in (|14[) . 

The subgraph protocol, on the other hand, depends on 
the size of the subgraphs. There are two extreme strate- 
gies, depicted in Fig. @] create subgraphs of degree 1, 
and connect them successively (SI), or create subgraphs 
of degree 2 every second node, and connect them at their 
common neighbor (S2). One could also adopt an inter- 
mediate strategy, in which the degree of each GHZ is 
selected at random, resulting in a mix of subgraphs of 
degree 1 and 2 — its performance, which will not be re- 
ported here, falls in between the two extreme strategies. 
In the first case (SI), Eq. ([22) reads 

(K x ) = 11 (1 p2 f{xu,x u+1 ,x u+2 ) (1 _ pi) *„©*„ +1 

UEV 

{{Kt) x -mi +1 r^ Xu+l )i (39) 

while in the second (S2) it is 



(K x )= 11 (i- Pa )«(»«.»«+i.*«+")(i-p 1 )*«+iffl*»+» 

((KtTHK*:_j x ^-HM :+1 r" +i ) . (40) 



Note that in the latter the product is over 

U = 2,4,- •• G Kven- 

Let us now define a domain of x = {x\x 2 . . . xjy) 
as a sequence {x u x u +\ . . .x u +i) where all x v , for v = 
u, . . . ,u + I, have the same value (either or 1), and 
where x u -\ and x u +i+\ have a different value. We can 
differentiate between a sequence of zeros or ones. Taking 
Eq. (131)1) as an example, one can observe that except for 
the term (1 — p 2 ), (K x ) factorizes into expected values 



in each component of ones, as the exponents of the noise 
terms can be different from zero if at least one of the bits 
x is equal to one, and the terms are uncorrelated with 
those on different domains. Considering all noise terms, 
Eq. (|39|) can be expressed as a function of the total num- 
ber of ones in the sequence x, n = \x\; the number c\ 
of domains of only one zero; and the number of domains 
C2* with two or more zeros: 



(41) 



(K x ) = (1 -p 2 ) Cl+2c -+»(l _ Pl )3<<*-H**) 

(KtT +C2 * (Kt) n - Ci - C2 * <ie^:) cH " 

The fidelity is then 

Fn = Jn 9{cuc 2 *,n,N) {K a .(c 1 ,c 2 *,n)) , 



where g(c\, c%* , n, N) is the number of (K x {c\,C2* ,n)) 
elements with these parameters in a graph of N vertices. 
This sum can be computed exactly by turning to generat- 
ing functions. The function generating g(c\, c-i* , n, N) is 
G(x, yi ,y 2 ,z) = J2g(c 1> c 2 *,n l N)x n y'iy?'z N . We de- 
rive its closed form in Appendix [B] which reads 



G = 



1 - xz 2 (l -yi + %i - y 2 )z) 
1 - (1 + x)z + x(l - yi)z 2 + x(yi - y 2 )z 3 

The fidelity for a state of size N is then 
1 1 



(42) 



(43) 



evaluated at 



x = (I - P2 ) (Kt) , (44) 
yi = (1 - P2 )(l -pi) 2 (K%) (K^)- 1 (K^KH) , (45) 
V2 = (1 - P2 ) 2 (l -Pi) 2 (K£) (Ktr 1 (K*>K%) , 

(46) 

and z = 0. The factorization of (K x ) for Subgraphs 
S2 and Bipartite protocols depend on slightly different 
parameters, and the corresponding generating functions 
can also be computed (see Appendix |B|) . 

We will first consider the ideal case of perfect purifi- 
cation, and later tackle the noisy case. In the ideal sce- 
nario, we can substitute all subgraph correlators by 1 and 
Eq. (T4"3l gives the exact fidelity for all p. The correspond- 
ing decay exponent /jv is plotted in Fig. [5] 

We now compare these exact results with the mean- 
field approximation. As we introduced in the previous 
section, (K x ) ~ exp[— /3H(x)]. For the Subgraph SI 
protocol, H(x) = J2uev h u{x), with 



h u {x) = d(x u , x u+ i,x u+2 ) + x^ 



(47) 



Recall that we are still considering perfect purification. 
Substituting x u — > s+S u and keeping only linear terms in 
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Figure 5: Rescaled decay rate /jv of a closed linear cluster 
of N = 10 nodes with perfect purification. Blue (squares): 
Bipartite A; red (diamonds): Bipartite B; black: Subgraphs 
SI (circles) and S2 (triangles). Dots correspond to the mean- 
field approximation, lines to the exact result using generating 
functions. 



S u , this takes the form of h^ F (x) = a u + b u x u + c u x u+ i + 
Xu+2, with 



(5-2s)s 2 
2 - 4s + s 2 
2 - 4s + s 2 

2 



d u = (1 - s) 



(48) 
(49) 
(50) 
(51) 



All qubits u are equivalent, so the coefficients a, b, c, d 
are the same for every node and the total Hamiltonian is 
h mf = AN + BJ2 U x u , with A = a and B = b + c + d. 
Hence, the fidelity in the mean-field approximation is 



7 MF 
N 



= e -PA" e -PBN/2l osh PB 



N 



(52) 



and the decay rate 



MF 



-log cosh- 



er 1 - (53) 



The same approximation can be made for the other 
protocols. The local Hamiltonian of Subgraphs S2 is 



h u (x) = 9(x u ,x u+ i,x u+ 2) +x u+ i ®x u+2 . 



(54) 



Note that, in this case, the total Hamiltonian is the sum 
for u even. In the case of Bipartite A and B, the local 
Hamiltonians are 

h u (x) = 2x u + x u ^i © x u+ i + 29(x u -i,x u ,x u+ i) (55) 
and 

h u {x) = x u + 6(x u -i,x u ,x u+1 ,x u+2 ), (56) 



respectively. The mean-field results are shown in Fig. [5] 
together with the exact result. The agreement is re- 
markably good, specially considering that this is a one- 
dimensional network configuration. We observe that the 
decay rate has a linear dependence in p but soon higher 
order (non-linear) terms start to kick in. Here, we have 
studied and plotted the solution for a wide range of p. 
This has been done for completeness and to check the 
validity of the mean-field approximation, but recall that 
the aim in this paper is solely to compute the decay rate 
to leading order in p. The reason being that in realistic 
scenarios p will be strongly limited by threshold values 
required for subgraph purification. 

We now consider the realistic scenario following the 
same procedure as before but taking into account the 
corrections due the noisy purified subgraphs. This can 
be done by approximating the correlators by their lin- 
ear correction around unity (|14|) . Consequently, an addi- 
tional term (f24| must be added to the local Hamiltonians 
g7j), ([54]) . ([55)1 . (15B1 for each of the different protocols. 
It is important to bear in mind that all subsequent re- 
sults have to be taken consistently up to leading order in 
p. This leading order is the one that enters in the fea- 
tures that we seek: the linear dependence of the fidelity 
decay rate around p ~ (for arbitrarily large N). From 
equation (f34|) and (|35|) it is immediate to see that at this 
order the whole contribution to the decay rate is given 
by the first two terms in evaluated at s = 1/2, i.e., 



TN 



i 



s=l/2 



N2 



N 



(57) 



The second equality is straightforward to get from the 
definitions of A and B, and it states that in this regime 
(low p) the decay rate is dominated by the exponents 
of the typical sequences x. The expected value of H(x) 
(over sequences x) can be easily carried out since x u are 
independent variables of mean 1/2 and it leads to fideli- 
ties that decay as p/nN with Jn equal to 9/2 = 4.5, 
43/16 w 2.7, 21/8 w 2.6, and 13/8 w 1.6 for the Bipar- 
tite A and B and Subgraphs SI and S2, respectively. 

These coincide with the decay rates obtained by the 
generating function method in the same limit, where in 
Eq. dm-flm) we substitute = (K^K^) =l-2p 

and {K^) C1+C2 ' = 1 - p. The results (see Fig. [BJ 
show that the subgraph protocols provide output fideli- 
ties comparable, if not better, than those given by the 
protocols based on channel purification. This is remark- 
able keeping in mind that purifying channels is much 
more demanding in terms of resources and efficiency. 

As we mentioned at the beginning of the section, the 
order of the CPHASEs in the Bipartite protocols can give 
different results, but the decay rates do not change much. 
For example, if the CHPASEs are first applied to every 
second edge, and then to the remaining edges, the decay 
rate is 143/32 w 4.5 for Bipartite A, while it does not 
change for Bipartite B. 
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Figure 6: Fidelity decay rate p/jv of a closed linear cluster 
of N = 100 nodes. From upper to lower lines: Bipartite A 
(blue), Bipartite B (red), Subgraph SI (black) and Subgraph 
S2 (black). Solid lines correspond to the generating function 
result, using the first order approximation of the purification 
scheme of Eq. (|14|l . Dashed lines are the first-order result of 
Eq. (I5T[) . Dots correspond to the mean-field approximation, 
dashed lines to the fidelity at p — > 0. 



B. Complex networks 

Now we study the behavior of the protocols in com- 
plex networks. Complex networks are characterized by 
statistical properties, and can be modeled as an ensem- 
ble of graphs Q, with a probability P(G) assigned to ev- 
ery graph G in the ensemble. A property O of a com- 
plex network is defined as its average over the ensem- 
ble, O = X^Gee P{G)0 G (note that the overline here 
stands for the ensemble average) . Some of these proper- 
ties are self-averaging, meaning that for large systems (in 
the limit N — > oo) a property of a given graph realization 
G is the same as the average over different realizations of 
G eG H3,|i||. As we discussed in Section lfflDl the free 
energy is typically self-averaging, and is the one we use 
here as a figure of merit. 

We consider random networks with uncorrelated de- 
gree distribution p k , which defines the ensemble. This 
family of networks includes the Erdos-Renyi model (when 
Pk is a Poisson distribution) as well as exponential and 
scale- free networks 32] . It is useful to express the degree 
distribution by its generating function g p (z) — '^2 k PkZ k . 
Related to p k is the probability that, following a random 
edge, we arrive at a vertex with k other edges. This prob- 
ability is = (k + l)pk+i/ (fc), and fc is called the excess 
degree. It is generated by g r (x) = g' p {x)/ (fc). The edges 
in the network are undirected, but the creation of the 
graph state via the Subgraph protocol is "directed-like" , 
as each node can have incoming and outgoing neighbors. 
Thus, we also consider the degree distribution pij, where 
i is the in-degree and j the out-degree. This distribution 
defines the implementation of the protocol. For exam- 
ple, in the linear cluster, SI was defined by pi.i = 1 



for all nodes, and S2 by p ,2 = 1 and p2,o = 1 at even 
and odd nodes, respectively. The distribution pij is 
constrained by J2i jti ~ i)Pi,j = 0; which means that 
(i) = (?) = (k) /2j but note that, in general, it can be 
correlated (pi.j ^ PiPj)- The function that generates pij 
is g p {x,y) = V, , />,.,.<•'//'. 

In the case of complex networks, we are interested 
in the average /at = J^ceg P(G)f G over the ensemble. 
We expect that, for large N, Jn goes to / due to self- 
averaging. In order to average /jv we have to compute 

' Geg x 

which corresponds to a quenched average, where the dis- 
order corresponding to the network topology, given by 
P(G), is frozen with respect to that of the correlation 
operator index, x. Computing this quantity is in gen- 
eral an extremely challenging problem, specially in our 
case that we have structured disorder and k-body inter- 
action between the spins. Nevertheless, we are interested 
in the particular regime of high temperatures (low p) far 
away from critical phenomena, long range correlations, 
and other difficulties that appear at low temperatures. 
In this regime we can directly use (|57l) : 



Geg x x 

(58) 

The Hamiltonians for the Bipartite protocols are of the 
form 

H G (x) = K{x) + h uAx). (59) 

uev (u,v)ee 

In the Bipartite A, the local Hamiltonians are 



h u {x) — x u + x v ■ x u + 2 x v 

+ 9(x u ,xtf u ) + x u + (J) x v (60) 

and 

while in the Bipartite B, 

h u {x^ — x u x v + 2x u + x u (62) 

and 

^W,tj(«^) — ^ | X u ,X Vl X w , X w ,Xj^^Xj^^ J . 

(63) 



11 



As before, the tilde over the summatory and the neigh- 
borhood stands for the order in which CPHASEs are per- 
formed in the local graph. We can substitute it by the 
expected effect of an edge, which depends on the number 
of edges already connected to nodes u and v. This num- 
ber is n u + n v with probability 2k J+k v Q£) (^)- Here, 
k u and k v are the excess degrees of the vertices in edge 
(it, v), so the average of the network ensemble has to be 
performed using probabilities r^ u and r^ v : 



h u ,v( x ) = ^n^2^2 r k u r k v KA x )- 

x k u fe„ 



(64) 



In the Bipartite A, the average effect of each of these 
edges is 1 — i [g r (3/4)] 2 . In the Bipartite B, it is 1 — 

[g r (l/2) + <? r (3/4)] 2 . The terms in h u (x), on the other 
hand, depend directly on the degree, and the average is 
performed over pk- 



K(x) 



2 N 



(65) 



Considering all these terms, and that the summation 
over V contains N elements, while that over E contains 
(k) N/2, the decay rates are 



fN = 



15 



1 



;9p(0) ~ o<? P (l/2) 



(fc) 
2 



l-\[g r {i/A)f 



(66) 



for Bipartite A and 



fN = \~ \g P (0) + & (l- 1 [ 3r (3/4) + . 9r (l/4)] 2 ) , 

(67) 
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Figure 7: Fidelity decay rate p/jv in an Erdos-Renyi network 
of mean degree (k) = 2 (top) and 4 (bottom). Bipartite 
A (blue, circles), Bipartite B (red, squares), and Subgraphs 
(black, diamonds). Bipartite A (upper, blue), Subgraphs 
(middle, black) and Bipartite B (lower, red). Dots are the av- 
erage over 10 random realizations of networks with N = 100, 
randomly sampling 1000 configurations of x. Dashed lines 
are the results at first order in p, Eqs. (|66|l . (|67|l . and (1681) . 



for Bipartite B. 

In the Subgraph protocol, where H G (x) 
Smgv h u (x), the local Hamiltonian is 



h u (x) 



+Xu(j (u> + 1) (1 -S 0jM ) + x u S Q>i ( U )S 0ij (u) 

+ S ' [^(^C u , X^fout , X Vl , . . . , X Va ) + X u © %Va\ • 

(68) 

Here Af^ is defined as in Eq. (f22j) . In this case, to average 
the Hamiltonian we have to take into account the directed 
degree probability p ij} h(x) = Y<ijPi3 h ( x )i g ivm g 



These results are valid for any network with uncorre- 
lated degree distribution pk- To compare the behavior of 
the three protocols we consider the Erdos-Renyi model 
(49Tj5l| , which is the maximally random graph under the 
only constrain that the mean degree (k) is fixed. The 
degrees in this model follow a Poisson distribution, so 
g p (x) = g r (x) = e' fc '^ x_1 ^, and for (k) > 1 there is a giant 
connected component of size comparable to the size of the 
network. The degree distribution pk is fixed by the net- 
work model, but in the Subgraph protocol one can play 
with the distribution pij (as long as J2i j Pi.j^i+j.k = Pk)- 
For simplicity, here we consider that the direction of each 
edge is selected at random, 



%3 



Pi+j (i + 3 



(70) 



^=! + ^ + jfc(0,0)- f<? P (l, 0) 
1 



1/2) -^(1/2,1/2)]. (69) 



and hence g p (x, y) = e^H - ^ -1 ). Fig.[7]plots pJn for the 
three protocols, showing that the best protocol depends 
on the mean degree of the network. It also shows an 
average over 10 random realizations of an Erdos-Renyi 
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graph of TV = 100, approximating the fidelity by the 
average of a random sample of 1 000 configurations of 
F ~ Mra Ex aamplo <^ samplo ). This approximation is 
valid in the low p regime, where the fidelity is dominated 
by the typical values of x. We observe that, in the Erdos- 
Renyi, the Subgraph protocol with random edge direction 
give a better (lower) decay rate /jv for (k) < 2.8, while 
above that the Bipartite B beats it. 

V. CONCLUSIONS 

We have proposed a protocol to generate a graph state 
that spans a distributed network of any topology, in the 
presence of noise in the communication channels and in 
the local quantum operations performed in the network 
nodes. The protocol distributes and purifies small sub- 
graphs, that are then merged to mimic the network struc- 
ture. This allows the use of multipartite purification for 
small systems (that depend on the degree of the nodes, 
rather than on the size of the network), which makes the 
protocol scale efficiently. We have compared this proto- 
col, which uses multipartite state purification and graph 
merging, to two other protocols that rely on channel pu- 
rification (quantum repeaters). 

We have benchmarkcd the protocols using as a figure of 
merit the fidelity of the generated graph state, or, more 
precisely, its decay rate, as the size of the network in- 
creases. Using generating function methods we have been 
able to compute exactly the fidelity for linear clusters of 
arbitrary size, allowing for a direct comparison with the 
approximate methods we develop. We have rephrased 
the problem of computing the effects of noise in the op- 
erations in terms of the thermal properties of a classical 
spin system, with the same interaction patterns as the 
underlying graph. Indeed, the fidelity itself can be seen 
as the analog of the partition function of such system, 
while its decay rate plays the role of the free energy. The 
well-known methods from statistical physics, such as the 
mean-field approximation, can be used to study its be- 
havior. 

We have also studied the three protocols in networks 
with a complex structure. These complex networks are 
present in many real-world systems, and are of particu- 
lar importance in the description of communication net- 
works. The subgraph protocol is specially motivated to 
be implementable in networks of this kind, as it is highly 



adaptable to any network topology and it only requires 
the nodes to have local knowledge of the network. 

Our results shows that the protocol using subgraph pu- 
rification and merging is comparable, and in some cases 
even better, than those which rely on bipartite states. 
This is quite remarkable, as the latter rely on quantum 
repeaters and require much more resources. In the com- 
plex networks case, the multipartite protocol could be 
possibly enhanced by devising an optimized strategy via 
the directed degree distribution pij. As an example, in 
networks with many nodes of degree 1, one should go for 
a directed degree distribution in which those nodes are as 
much as possible the receivers of GHZ states, instead of 
the senders, so less connections would need to be made. 
In other words, one should tend to pi t o as high as possible 
(compared to po.i)- Also, in nodes of too high degree, the 
protocol might fail due to the noise threshold in the mul- 
tipartite purification, which depends on the size of the 
GHZs (and thus, on the degree of the network nodes). In 
this case, one could always separate the node in two (or 
more) and treat them as independent nodes, each creat- 
ing and distributing a subgraph among a subset of the 
neighborhood. This might prove useful in networks with 
a scale-free degree distribution, which have a long tail 
and a high presence of hubs. 

It remains an open question to relate the decay in fi- 
delity with the actual use one can make of graph states. 
Clearly, there is a regime where the fidelity is exponen- 
tially small, but a finite decay rate still signatures valu- 
able resources. It is important to emphasize that for 
the three protocols, as well as for all the networks un- 
der study, we have not only computed this decay rate, 
but also completely characterized the noise graph states 
that emerge form the protocol. For a particular applica- 
tion, having full knowledge of the generated graph gives 
probably a closer picture of the generated resource. 
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Appendix A: Multipartite purification protocol 

The recursive protocol 38] is an entanglement purifica- 
tion protocol for multipartite, two-colorable graph states, 
which has also been extended to general graph states [3!| . 
It operates on two identical copies 

Xl,X 2 

of a two-colorable graph state of size N (with colors A 
and B), and consists of two subprotocols (PI and P2), 

each of which purifies one of the two colors. Here K&} is 
a stabilizer element of state pi. In each subprotocol, in- 
formation about p\ is transferred to p2- Then, p\ is kept 
or discarded depending on the outcomes of measurements 
on p2- In PI, a CNOT2^i is applied to every node in A, 
and a CNOTi^2 to every node in B. This transforms the 
state to 

— V I I K [Al) K {Bl) K (B2) 

x\ ,x 2 

where a and b are the elements of x that correspond to 
colors A and B. The stabilizer K^^j a2 corresponds to 
qubits of pi in color A, and the modulo 2 summation of 
the index a\ © a-i is made elementwise. To keep notation 
short, we contract K^^K^ = K^^ ^, where 
the first subindex corresponds to color A and the second 
to color B. Then, every node in p 2 is measured: X in 
nodes of color A (outcomes £) and Z in nodes of color B 
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(outcomes £). This gives the (unnormalized) state 



E(- 1 ) (ceec) - a2 « ) , bl >« ) , bl X ) 



iffia2,bi 



where £ © C stands for © ® v& j^ Cu fo r all u G A 
The state is selected if £®® C = 0. Summing over all the 
possible outcomes, the final (post-selected) state after PI 
is 



7\A\ 1^ \ K a 1 K b 1 j^a^b, ) 2 K ai ® a2 K bl 



2 N \ 



P2 is equivalent, with colors A and B interchanged. 

For our protocol, we need the fixed point for a GHZ of 
size j + 1, with a central node colored as A and j leaves 
colored as B (\A\ = 1 and \B\ — j). Noise can come from 
CNOTs as 1 — p2 and from measurements in state pi as 
1-Pi: 



a 2 =0 

(1 -pjj^^es 6 ) JJ(1 _ p2 )«(i^2,b) 

<^> (P2) =4 E (*S X, *&.) 

6 2 =0 

(1 _ p2 )«(«.©6 6 s2 , 2,© il6S 6) J|(l -p 2 ) e ( a ^2,b) 

6( 



6G-B 
\6a 



fees 



We now consider pi = p 2 = p and approximate (K a K b ) 
at first order in p. Let the unnormalized {KaK^y 1 ' ~ 



1 - /8 ,|6|P and (K a K b 



(P2) 



1 — "a,|6|P- Composing PI 



and P2 we can find the fixed point at first order in p. In 
PI, each P a ,\b\ equals So, 161 +51,161 plus a constant term: 



h),\b\ = «0,|6| + ai,|b| + j + 1 



5 l,|b| = «0,|6| + "1,161 + 2^ + 1 )- 



These terms are normalized dividing them by (1) = 
(KqKq), so the normalized first-order coefficients for PI 
read @ a M = P a ,\b\ ~ Po,o- Similarly, for P2: 



"o,|6| 



E^ ^ + Pom) + J + 1 + 



After normalization (dividing by (KqKq)), 



(K K b ) ~ l 



(XiKb) ~ 1 - (j + l)p. 



The fidelity is 
F~ 1- 



2-'- 



tE 



b=0 



+ U + 1) 



= i--(i + i)p. 



Appendix B: Generating functions for the domains 
in the closed linear cluster state 

In this appendix we derive the generating functions 
from which the fidelity of a closed linear cluster can be 
derived using Eq. (14U)) . In all cases, errors from CPHASEs, 
measurements and multipartite purification are tracked 
by error parameters p 2 , p\ and p, respectively. 



1. Subgraphs SI: g(n,ci,C2*,N) 

In a closed linear chain, nodes can have values "0" or 
"1" . In this chain, a domain is a sequence of adjacent 
nodes with the same value (0 or 1), surrounded by nodes 
of different value. In a chain of N nodes, let n be the 
total number of ones, c\ the number of domains of ones 
preceded by a domain of only a zero and c 2 * the number 
of domains of ones preceded by a domain of at least two 
zeros ( C\ + c 2 * is the total number of domains of ones). 
Recall Eq. (|4"Tj) , which gives the correlator of the protocol 
Subgraph SI in terms of n, c\ and c 2 *\ 

(K m ) = (1 -p 2 ) Cl+2c -+"(l - Pl f^+^) 

Then, g(n, ci, c 2 *, N) is the number of different configu- 
rations of that chain with given parameters, and 



G(x,y 1 ,y 2 ,z) 



g(n,c 1 ,c 2 ,N)x n y?y?*z N 



its generating function. Each variable x, y\, y 2 , and z 
"counts" the number of ones, domains of one zero, do- 
mains of two or more zeros and the total size, respec- 
tively. The function G can be found by joining simpler 
distributions. We can think of the linear chain as a con- 
struction of domains of zeros and ones joined together. 
Consider the set of domains of zeros, {0,00,000,...}, 
each domain of a given size appearing only once. The 
number of domains of size N in this set is (In = 1, which 
is generated by 



D(z) 



n>l 



dftZ 



N 
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We can differentiate between the set of domains of only a 
zero, {0}, and that of two or more zeros, {00, 000, . . . }. In 
this case, the generating functions are respectively z and 
D{z) — z. The set of domains of ones, {1, 11, 111, . . . }, is 
generated by the same function, but here each element 
contributes to the total size of the chain and to the num- 
ber of ones. Its generating function is thus D(xz). The 
function generating the set of pairs of domains, the first 
of zeros and the second of ones, is 

P = P(x,y 1} y 2 ,z) = {yxz + y 2 [D(z) - z]}D(xz). 

The function for all possible combinations of pairs of do- 
mains of zeros and ones (in order) is 

— - = 1 + { Vl z + y 2 [D{z) - z}} Y^-£. 

Here, the first element (1) counts the case where there is 
no pairs at all, and the second, Pj (1 — P), to that where 
there is at least one pair. We can now add nothing at all, 
a possible domain of zeros at the end, a possible domain 
of ones at the beginning, or both. The final generating 
function is then 



G= \l + [ yi z + y 2 (D(z)- : 

+ l + y 2 D{z)D(xz) 



)]£M}[1 + D{xz)+P\ 

1 



Note that when we added a domain of zeros (the term 
with D(z)), we changed the variable j/i for y 2 , to take 
into account that the domain of zeros is now of size 
greater than one (because we are considering a closed 
linear chain). Simplifying, we obtain 



G 



1 - xz 2 {\ ~yi + 2(?/i - y 2 )z) 
l-(l + x)z + x(l - y^z 2 + x(yi - y 2 )z 3 ' 



(Bl) 



D(z). 



1. Bipartite B: g(n, ci, C2, C3* , N) 

In Bipartite B protocol, the correlator depends on sim- 
ilar parameters, but here we have to differentiate between 
the number of domains of zeros with one element (ci), 
two elements (02) and three or more elements (03*): 

(K x ) = (1 -pi)"(l _p 2 ) C i+2c 2 +3c 3 .+n 

(l-2p)"(l-p) 2c2+2c3 *. (B2) 

This can be achieved by a small modification of the pre- 
vious generating function. Now, 

P = P(x,y 1 ,y 2 ,y 3 ,z) 

= { Vl z + y 2 z 2 + y 3 [D(z) - z - z 2 ]} D{xz), 



and 



G 



(1 + xz 2 {-\ + yi - 2y lZ + 2y 2 z + 3(-y 2 + V^)z 2 )) 
(1 + z{-l + x{-l + z(l + 2/1 (-1 + z) + z{y 2 (-l + z)- 2/3-2))))) ' 

I 



(B3) 



3. Bipartite A: g(n, ci, C2» , C2* , N) 

In the Bipartite A, we also need to count the number 
c 2 * of domains of two or more ones. The correlator in 
this case reads 

(K x ) = (1 _p 1 )2n+2c 2 .+2c 2 .( 1 _ p2) 2c 1+ Ac 2 +2n 

(1 - 2p) 2c2 '+ 252 * (1 - p) n - 2S 2* . (B4) 



ated by xz and D(xz) — xz, respectively. The extended 
function generating the set of pairs of domains is now 



P = P(x, yi,y 2 ,w 2 , z) = {y 1 z + y 2 [D(z) - z]} 

{xz + w 2 [D(xz) — xz]} . 



We count c 2 * using variable w 2 , and differentiating be- 
tween the sets {1} and {11,111,...}, which are gener- Proceeding as in the previous case, we obtain 



G 



1 + xz 2 {-\ + yi- 2 Vl z + 2y 2 z) + (-1 + w 2 )x 2 z 3 { yi (2 - 3z) + 3y 2 z) 
1 + z{-\ + x(-l + z(l + (1 + (-1 + w 2 )xz){ yi {-l + z)- y 2 z)))) 

I 



(B5) 



4. Subgraphs S2: 3(7101, mo, mi, ci, c r , N) 



a; as a sequence of elements 00, 01, 10 and 11 (the first 



In the Subgraphs S2 protocol, the sum is performed 
over even nodes. In this case, it is convenient to express 
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digit corresponding to an odd node, and the second to an 
even node). Each 01, 10 and 11 contribute to one CPHASE 
noise, as well as each domain of 00 which is preceded by 
a 01 or a 11. Moreover, each 01 and 10 contribute to one 
Y measurement noise. Finally, each 01 and 11 contribute 
to a 1 — 3p noise of the purified subgraph, each 10 to a 
1 — p and each domain of 00 followed by a 10 or a 11 also 
to a 1 — p. Thus, we need the number of configurations 
with noi, nio, and tin number of 01, 10 and 11 elements, 
and c ; and c r domains of 00 preceded by 01 or 11 and 
followed by 10 or 11, respectively. The correlator is 

(K x ) = (1 -pj^oi+nio+mi+c,^ _ pi )n 01 +n 10 

(1 - 3p) noi+ni1 (1 - p) ni0+Cr . (B6) 

Now, each element contributes with z 2 to the size of 
the chain. The domains {00, 00 00, ... } are generated by 

2 

x z _ 2 . A domain made of elements 01, 10 and 11 of any 
size (including 0) is generated by 



1 



1 - (a:oi + £io + x\\)z 



2 • 



and one which ends (or begins) with, say, element 01 (and 
thus is of size at least 2) is generated by 



i 7 i i r^ x oiz . 

Proceeding as in the previous cases, the generating func- 
tion is 



G 



(a?oi +x w + x n )z 2 z 2 
1 - (a:oi + x w + xn)z 2 1 - z 2 



G 



1 



where P is 



1-P' 
(B7) 



{ {xoiyi + x 10 y r + xnyiy r )z 2 



1 - (xoi + xio + xn)z 2 

[{x w + x n )y r + xoi}} - - 2 
i — z 



z [(x i + xu)yi + xio] 

^2 



